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ABSTRACT 

We present VLT intermediate resolution spectroscopy of UY Vol, the optical counterpart 
of the LMXB X-ray burster EXO 0748-676. By using Doppler tomography we detect narrow 
components within the broad Hell A4542, A4686 and A5412 emission lines. The phase, ve- 
locity and narrowness of these lines are consistent with their arising from the irradiated hemi- 
sphere of the donor star, as has been observed in a number of LMXBs. Under this assumption 
we provide the first dynamical constraints on the stellar masses in this system. In particular, we 
measure K 2 > K em = 300 ± 10 km s -1 . Using this value we derive 1 M Q < M 1 < 2.4 M 
and 0.11 ^ q ^ 0.28. We find Mi ^1.5 M© for the case of a main sequence companion star. 
O ur results are consistent with the presence of a massive neutron star as has been suggested 
bv lOzell d2006[) . although we cannot discard the canonical value of ~ 1.4 M Q . 
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1 INTRODUCTION 

EXO 0748-676 (=UY Vol) is a Low Mass X-ray Bi nary (LMXB) 
disco vered as a transient source by EXOSAT in 1985 jParmar et al.1 
1985). Instead of returning to the quiescent state the system had 
remained in an active state ever since, showing typical properties 
of persis tent LMXBs. Howeve r, it has returned to the off-state very 
recently drlvnes & Jonesll2008l) . 

EXO 0748-676 also exhibits irregular X-ray dips and periodic 
eclipses which have provided an excellent determination of the or- 
bital period (P=3.8 2 h) and the orbital inc lination (i = 75° - 82°; 
IParmar et all 1 19861) . On the other hand, IVillarreal & Strohmaveil 
1 2004) have detected a burst oscillation at 45 Hz during type I X- 
ray burst events which is associat ed with the neutron star (NS) spin 
frequency. Recently, Ozel (2006) have suggested the presence of a 
massive NS with Mi = 2.10 ± 0.28 M and i?i = 13.8 ±1.8 km 
in EXO 0748-676. This result is based on the gravitational redshift 
measur ed in O and Fe abso rption lines arising from the surface of 
the NS l lCottam et al.l20 02) and the assumption that the strong ther- 
monu clear X-ray bursts displayed by the system dGottwald et al.1 
Il986h can be modelled by a phase of symmetric expansion of 
the NS radius and a cooling ph ase where the X-r ay emission is 
Eddington-limited. In contrast, IPearson et alj fc006l) have reported 
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He II A4686 Doppler tomograms which show a spot consistent with 
the expected gas stream position. These authors favor a NS of 
1.35 Mq by fitting different gas stream trajectories to the position 
of the spot. 

In this paper we present VLT-FORS1 spectroscopy of the 
V ~ 17 optical counterpart (UY Vol) of EXO 0748-676 which 
reveals the presence of narrow emission components within the 
He II emission lines. In recent years, narrow emissi on lines from the 
donor have been det ected in severa l LMXBs (see lCornelisse et ail 
l2008lfor details, alsolCasares|2009l) since they were first discovered 
by ISteeghs & Casares! J2002I) in Sco X- 1 . The Roche lobe filling 
donor star intercepts the energetic photons from the inner accre- 
tion disc resulting in the observed optical emission lines from its 
surface. The Doppler motion of these lines traces the orbit of the 
companion star, thereby providing the first constraints on the bi- 
nary parameters for persistent systems where the photospheric light 
from the companion is otherwise swamped by the bright accretion 
disc. It is remarkable that this technique has prove d the presence of 
a black hole with a mass > 6 M Q in GX339-4 favnes et al.1l2003l : 
Munoz-Darias et al. 2008), presented evidence for a massive NS 
in Arjl X-l jCornelisse etaUfcOOTl) and X1822-371 JCasares et ail 
l2003l :lMunoz-Da rias et al.ll2005l) . and provided the first contraints 
on the binary parameter s in six other LMXBs, inc luding the proto- 
typical system Sco X-l dSteeghs & Casa res 2002). In this letter we 
present the first dynamical constraints on the stellar masses in EXO 
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Figure 1. Average spectrum of EXO 0748-676 after combining 21 expo- 
sures obtained with good seeing conditions (< 0.6 arc-sec) and out of 
eclipse. The most prominent lines are indicated for clarity. 
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Figure 2. Diagnostic diagram of He II A4686 obtained by fitting eq.[T]to the 
radial velocity curves extracted for each value of the Gaussian separation a. 
The vertical dashed line shows the selected value for a and the horizontal 
lines the values adopted for each parameter. 



0748-676 based on the study of the narrow components present 
within the He II A4541, 4686 and 5412 emission profiles. In section 
2 we describe the observations and data reduction process, whereas 
the data analysis is presented in section 3. Results are discussed in 
sections 4 and 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed EXO 0748-676 on 22 Jan 2008 using the FOcal 
Reducer and low dispersion Spectrograph- 1 (FORS1) attached to 
the VLT-UT2 (Kueyen) telescope at the Observatorio Astronomico 
Cerro Paranal (Chile). We used the GRIS_1200g holographic grat- 
ing obtaining a total of 28 spectra of 600s each. The slit width was 
fixed at 0.51-arsec resulting in a spectral resolution of ~ 90 km s~ 
and a wavelength coverage of AA4185 — 5600 A. The seeing stayed 
between 0.4-0.6-arsec during the first 4 hours but it increased up to 
1-arsec during the last hour of the run, making slit losses significant 
and absolute flux calibration unreliable for this part of the night. In 
total we covered ~ 1.3 binary orbits (i. e. ~ 5 h of data). The im- 
ages were de-biased and flat-fielded and the spectra extracted using 
conventional IRAF routines. The wavelength calibration was per- 
formed by fitting a fourth-order polynomial to an arc image taken 
during daytime. This resulted in a dispersion of 0.73 A pix -1 and 
rms scatter ^ 0.04 A. We used the sky spectra to correct for veloc- 
ity drifts due to instrumental flexure which were found to be ^ 10 
km s" 1 (i. e. ^ 0.25 pix). The zero point of the wavelength scale 
was established from the position of the strong OI A5577.338 sky- 
line. The spectra were finally flux calibrated using observations of 
the flux standard LTT 1788 and exported into our analysis software 
(MOLLY). 
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Figure 3. Trailed spectra showing the orbital evolution of the emission lines 
He II A4542, A4686 and A5412 (from left to right) in 20 phase bins. Two 
orbital periods have been plotted for clarity. 



In this paper we will focus on the analysis of the He II emission 
lines where narrow components are found. In order to compute 
the orbital phase o f each spect r um w e used the constant-period 
ephemeris given by Wolff et al. (2002j). The zero-phase has been 
refined by fitting the mid-eclipse time of an eclipse observed by 
RXTE on 30 Jan 2008 (i.e. only 8 days after our observations). We 
find an error of ~ 0.012 in orbital phase which has been corrected 
in the ephemeris. We finally adopted: 



3 DATA ANALYSIS 

Fig. [T] shows the average spectrum of EXO 0748-676. In order to 
produce this figure we used the 21 out-of-eclipse spectra taken 
with good seeing conditions, comparable to the slit-width. The 
H7 and H/3 double peak emission lines together with the strong 
Hell A4542, A4686, A5412 features and the Bowen blend at 
AA4630 - 50 are clearly seen in the spectrum. He I A4922, A5016, 
and Fe II A5284 are also detected, although they are less prominent. 



T (HJD) = 2446111.57285(3) +0.1593378191(1) x N 
3.1 Systemic velocity 

As in our previ ous work, we decided to a pply the double gaus- 
sian technique JSchneider & Yound Il980|) to the wings of the 
He II A4686 emission line in order to obtain a first estimate of the 
orbital parameters of the NS. We have used relative Gaussian sep- 
arations between a = 800 — 2000 km s~ x in steps of 25 km s _1 . 



Dynamical constraints on the neutron star mass in EXO 0748-676 3 




-500 500 -500 500 -500 500 

Velocity X (km/s) Velocity X (km/s) Velocity X (km/s) 



Figure 4. Doppler maps corresponding to Hell A4542, A4686, A5412 (from left to right) computed using our refined ephemeris. A systemic velocity of 
7 = 70 km s — 1 was used for the three cases. We have over-plotted the position of the compact object, the gas stream trajectory and the Roche lobe of the 
companion star assuming K\ = 100 km s _1 and K2 = 360 km s — 1 . 



This range has been chosen to avoid low velocity components asso- 
ciated with the companion star and the outer disc. After fitting the 
expression 

V{4>) =7 + ^1 sin 2-K{<f> - <j> ) (1) 

to the radial velocity curves extracted for each value of a, we obtain 
the diagnostic diagram presented in Fig. [2] We adopt the value of 
a ~ 1700 kms -1 where the phase is correct and also happend 
to occur near the plateau in the diagnostic parameter a(K\) / K\ 
dShafteretal.1 19861) . This results in 7 = 30±10 km s" 1 and Ki = 
100 ± 20 kms" 1 . 

This method assumes that the emission line gas orbits the 
compact object and correctly traces its orbital dynamics. Non- 
isotropic emisison and/or significant disk asymmetries are known 
to lead to significant distortions in the radial velocity curves and 
thus the inferred K\, For this reason the values provided by the 
diagnostic-diagram should be taken as a rough approximation to 
the real ones, as will be discussed in the next section. 

3.2 Doppler mapping 

Extracting radial velocities of faint short-period systems is not 
straightforward. The situation is even more difficult when dealing 
with narrow spectral lines which require high spectral resolution. 
The spectral resolution also degrades with exposure time due to or- 
bital smearing, whereas the signal-to-noise ratio (S/N) obviously 
increases for longer integrations. Therefore, in order to avoid sig- 
nificant orbital smearing one has to use short exposures (600 s in 
our case) which results in data sets with good orbital sampling but 
poor S/N. 

Using Doppler tomography (Mars h & Hornell 19881) instead of 
standard radial velocity curves has the benefit of using all the spec- 
tra at the same time, allowing one to map the brightness distribu- 
tion of a weak emission line in velocity space. This method has 
been us ed in almost all prev ious studies of narrow emission lines 
(see e.g. lCasares et alj|2006l) . and only for Sco X-l (with V ~ 12 
and P=18.9 h) has the application of both techniques separately 
(i.e. radial velocity curves and Doppler mapping) been possible. In 
that case a good agreeme nt between the two methods was found 
dSteeghs & Casarej2002h . 

Fig. [3] shows the evolution of the He II A4542, A4686 and 
A5412 spectral lines along the orbit. An S-wave with a maximum 



radial velocity of ~ 500km s _1 at orbital phase ~ 0.5 appears 
in the 3 maps. This translates into bright spots at V(Vk,Vy)=(- 
500,200) km s _1 in the Doppler maps of the three He II lin es (Fig. 
[4^. Similar spots were also found bv lPearson et alj d2006h in sev- 
eral emission lines (e.g Hell A4686). However, in this study we 
also detect in the three maps compact spots at the position of the 
donor (i.e. Vx = km s _1 ). The velocities associated with these 
spots are in the range Vy = 270 — 339 km s _1 and are shown in 
tabled] 

We have computed Doppler maps for a wide range of 7 val- 
ues bracketing 7 = 30 kms -1 , the value favored by the diag- 
nostic diagram. We swept between 7 = — 50 and 100 kms -1 in 
steps of 10 km s _1 . The effect of using a wrong value of 7 when 
computing the maps is to obtain out-of-focus spots, which become 
compact when approaching the correct value. By fitting 2D gaus- 
sians to the spots as a function of 7, we favor systemic veloc- 
ities in the range 20 < 7 ^ 80 km s^ 1 . We finally selected 
7 = 70 kms -1 because it produces significantly more compact 
spots in the three maps. In any case, the velocity measurements (ta- 
ble[T|l derived with 7 = 70 km s _1 are consistent within the errors 
with those obtained with 7 = 30 km s _1 . 

The K em velocity of a given irradiation line depends on the 
disc thickness and how transparent it is for the irradiating pho- 
tons. The three Hell lines that we study here have similar ion- 
ization potentials and, hence, one should expect them to show the 
same K cln value. This is consistent within the error bars with the 
results that we have obtained. Therefore, we decided to take the 
weighted mean value of the three lines as our best estimate of K cnl , 
which results in K eul = 310 ± 10 kms" 1 . Hence, we will use 
K cul — 300 km s _1 as 1-sigma lower limit to K2. 



4 REPROCESSED EMISSION FROM THE COMPANION 

Using our K em value derived from the He II lines and the re- 
fined orbital ephemeris it is possible to remove the orbital veloc- 
ity from each spectrum and compute the average spectrum in the 
rest frame of the companion. Fig|5]shows this spectrum, where nar- 
row components standing out in the broad He II A4542, A4686 and 
A5412 emission lines are clearly seen. As an independent test to 
the Doppler analysis, we have computed the spectrum in the rest 
frame of the donor sweeping a wide range of Kem values. We find 
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Table 1. K em (V x , Vy) measurements (km s" 1 ) in EXO 0748-676. 

A4542 He II A4686 He II A5412 

(28, 300) ±18 (7, 325) ±14 (-14, 290) ±20 



that, for a range of K cm values around ~ 310 ± 10 km s" 1 , the 
narrow components show show the highest intensities together with 
narrow widths consistent with the spectral resolution of ~ 90 km 
s . Taking into account the position of the spots in the Doppler 
maps (i.e. orbital phase and velocity), the narrowness of the com- 
ponents and the presence of similar high excitati on features in all 
the p ersistent LMXBs that we have studied (see Corneli sse et al.l 
2008) we associate these features with X-ray reprocessing from the 
irradiated hemisphere of the donor. 

It is also noteworthy that we have a tentative detection of the 
N III and C III components within the Bowen blend in the rest frame 
of the companion star (left panel in Fig. [5}. The weakness of these 
lines is probably due to the reduced number of spectra. We only 
cover one orbit and at least two orbits of data are usually req uired to 
clearly detect these transitions (e.g. see lCasares e t al. 2006). How- 
ever, the He II emission lines are very strong in this system com- 
pared to others previously studied and show narrow components 
that provide us with the opportunity to measure K crn . It is espe- 
cially interesting to note that EXO 0748-676 is the only system in 
the Bowen survey where Hell A4542 is as strong as the Bowen 
blend. Hell A4686 emiss ion from the donor was also reported by 
ISteeghs & Casaresl {2002) in Sco X-l, with a K em value similar to 
that obtained from the Bowen lines. However, in other bright sys- 
tems where high signal-to-noise studies have been performed (e.g. 
X1822-371 Icasares et alJ fcoOS) no He II emission from the donor 
was found. 



5 SYSTEM PARAMETERS 

EXO 0748-676 offers important advantages with respect to other 
persistent LMXBs. The binary inclination is a known funct ion of 
the mass ratio owing to the X-ray eclipses ( iHvnes et al]|2006h . For 
instance, i — 82° yields q = 0.12, whereas i — 75° implies q = 
0.4. Moreover, since X-ray eclipses and X-ray bursts from the NS 
surface are observed, the NS must be visible, and hence the disc 
opening angle a, which obscures part of the irradiated face of the 
donor, is also related to the inclination through a ^ 90° — i . 
Therefore for a given q, the inclination angle and an upper limit to 
q are set. This is very useful in order to compute the K-correction. 
Emission lines formed on the inner hemisphere of the donor star as 
a result of X-ray reproc essing only provide a lower limit (K em ) to 
the true K2 -velocity. In lMunoz-Darias et ail ( l2005h we model the 
deviation between the light-center of the reprocessed lines and the 
center of mass of a Roche lobe filling star in a persistent LMXB, 
including screening effects by a flared accretion disc (i.e. the K- 
correction). We find that this correction depends on q, a and weakly 
on the orbital inclination. We now apply the K-correction to EXO 
0748-676 as a function of q, including the above limits on q and i. 



5.1 The mass of the neutron star 

Fig. [6] shows the lower limit to the NS mass as a function of q 
obtained from K cul — 300 kms -1 (lower solid line). Assuming 
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Figure 5. Normalized average spectrum of EXO 0748-676 in the rest frame 
of the companion (using K CIn = 300 km s — 1 ). Note the narrow compo- 
nents standing out in He II A4542, A4686 and A5412. 



K\ ~ 100 kms" 1 (favored by the diagnost ic dia gram) we ob- 
tain q < 0.35, in agreement with Hvn es et al] fc006l) . On the other 
hand, if we use the maximum K-correction (a = 0°) and the up- 
per limit Kcm = 320 kms" 1 we obtain an upper limit to the 
NS m ass (upper solid line in Fig. iMever & Mever-Hofmeisteg 
dl982h suggest that, even in the absence of irradiation, accretion 
discs in LMXBs have a > 6° , which is also consistent with the re- 
sults bv lD' Avanzo et al. (2006) for Cen X-4 in quiescence. Apply- 
ing the the K-correction with a J> 6° and K crn — 320 km s" 1 we 
obtain a less conservative upper limit represented by the dotted- 
dashe d line in Fi g. [6] 

IPearson et all l2006t) favored q = 0.34 y Mi » 1.35 M by 
fitting gas stream trajectories to the extended spot observed in the 
He II ADoppler map. This pair of values, which is marked by the 
dotted line in Fig. [6] is not consistent with our results. Con versely, 
our results are consistent with those derived by Ozel (2006), which 
are represented by the grey region in Fig. [6] Nevertheless, our val- 
ues are also consistent with the presence of a canonical NS with 
~ 1.4 Mq in EXO 0748-676. However, for this case to work we 
need q ^ 0.24 

5.2 The nature of the donor star 

The parameter constraint for EXO 0748-676 are consistent with a 
main sequenc e star and clearly a g iant star cannot fit within the 3.8 
h Roche lobe. Faulkn er et al. Ul972h showed that the mean density 
(p) of a Roche lobe-filling companion star is determined solely by 
the binary period P: 

p^H3P- 2 gcm- 3 

where is the orbital period in hours. For EX O 0748-67 6, with 
P h = 3.82 hours, we obtain M 2 ~ 0.42 M dCoxll2000l) . which 
defines the right dashed line in Fig. [6] through Mi = M2 /q. 
Combining this relation with the lower, thick solid line we obtain 
Mi 1.5 Mq for the case of a main sequence donor. 
Alternatively, the donor star could have evolved before the onset of 
mass transfer. For this to happen within a Hubble time, the binary 
must have experienced an episode of unstable mass transfer with 
q > 1 in a thermal time scale. The final product of this scenario is a 
LMXB formed by a donor star with an evolv ed nucleus transfering 
material onto a NS. ISchenker & Kind J2002I) have computed com- 
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panion star masses versus the orbital period for the limiting case of 
a nearly all helium nuclei. Using their relation for the case of EXO 
0748-676 we obtain M 2 ;> 0.16 M which yields the left dashed 
line in Fig. [6] Altogether, the allowed region in the Ml — q space of 
parameters yields 1 M < Mi < 2.4 M and 0.11 < q < 0.28 
whereas a plausible a ^ 6° implies Mi ^ 2.0 M . It also in- 
teresting to consider a limiting case with q — 0.28, a — 6° and 
Kem = 320 kms -1 together with K\ ~ 100 kms -1 from the 
diagnostic diagram. For this case we obtain q 0.22 and therefore 
Mi ^ 1.3 Mq, favoring the main sequence scenario. However, 
we note that the unstable mass-transfer scen ario was proposed to 
explain the ob servations in XTE Jl 118+480 iHaswell et al.ir2002h 
and Cyg X-2 {Po dsiadlowski & Rappaport l2000h . More recently. 
iRodrfguez-Gil et alj d2009h show that the mass of the companion 
star in the cataclysmic variable H S 0218+3229 is also in e xcellent 
agreement with the predictions o f lSchenker & Kind J2002h . 



6 CONCLUSIONS 

We present the first dynamical determination of the NS mass in 
EXO 0748-676. Assuming that the narrow components detected 
within the He II A4542, A4686 and A5412 emission lines arise 
from the irradiated hemisphere of the donor we measure A' cm = 
300 ± 10 kms -1 which results in 1 M Q Mi sC 2.4 M Q and 
0.11 ^ q ^ 0.28. On the other hand, the mass of the donor star 
is constrained to 0.16 M Q ^ M% < 0.42 M Q through unstable 
mass-transfer/main sequence models. For the case of a main se- 
quence donor we find Mi 1.5 Mq. 

Our results are consistent with the Ozel ( 20 03) scenario an d 
rule out other system parameters proposed by Pearson et al. |{2Q03). 
However, our observations are also consistent with a canonical 
mass NS and more observations are required to further constrain 
the mass of the NS in EXO 0748-676. Recently, EXO 0748+676 
has been re ported to reach quiesce nce at R = 21.9 after ~22 years 
of activity JHvnes & J ones 2008). This provides an excellent op- 
portunity to measure the orbit of the companion star and determine 
the mass of the NS. 
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